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Weak approximations have been developed to calculate the ex- 
pectation value of functionals of stochastic differential equations, and 
P> ' various numerical discretization schemes (Euler, Milshtein) have been 

studied by many authors. We present a general framework based on 
semigroup expansions for the construction of higher-order discretiza- 
i^H ' tion schemes and analyze its rate of convergence. We also apply it to 

, approximate general Levy driven stochastic differential equations. 

1. Introduction. Weak approximation problems play an important role 
in the numerical calculation of E[f[Xt{x))\ where Xt{x) is the solution of 
the stochastic differential equation (SDE) 



O 



O 

o 



Xt{x)=x+ f %{Xs-{x))ds+ I V{Xs^{x))dBs 

n 1N -^0 Jo 

OO ■ (1-1) ^ 

+ / h{Xs-{x))dYs 
Jo 



with smooth coefficients Vo-.K^ ,V = {Vi, . . . ,Vd),h:'R^ 

^ , whose derivatives of any order (> 1) are bounded. Here Bt is a d-dimensional 

5— ( ' standard Brownian motion and Yf is an d-dimensional Levy process associ- 
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ated with the Levy triplet (6, 0, z^) satisfying the condition 




for 



We denote briefly by E[f{XT{x))] - E[f{X^'> (x))] = ©(l/n'") the above 



situation, and say that Xj^ is a mth-order discretization scheme for Xt 

(n) 

or that Xj. is an approximation scheme of order m. The Euler scheme 
is a first-order scheme, and has been studied by many researchers. Talay 
and Tubaro [19] show the first-order convergence of the Euler scheme and 
second-order convergence with the Romberg extrapolation for continuous 
diffusions. The fact that the convergence rate of the Euler scheme also holds 
for certain irregular functions / under a Hormander type condition has been 
proved by Bally and Talay [3] using Malliavin calculus. For the general 
Levy driven case, the Euler-Maruyama scheme was first studied in Protter 
and Talay [18], see also Jacod and Protter [8] and Jacod et al. [7] (for 
smooth /). The Ito-Taylor (weak- Taylor) high-order scheme is a natural 
extension of the Euler scheme although is hard to simulate due to the use 
of multiple stochastic integrals. A discussion on the Ito-Taylor scheme with 
the Romberg extrapolation can be found in Kloeden and Platen [9]. 

In the continuous diffusion case, some new discretization schemes (also 
called Kusuoka type schemes) which are of order m>2 without the Romberg 
extrapolation have been introduced by Kusuoka [11], Lyons and Victoir [13], 
Ninomiya and Victoir [16], Kusuoka, Ninomiya and Ninomiya [12], Ninomiya 
and Ninomiya [15] and Fujiwara [5, 6] (m = 6). The rate of convergence of 
these schemes is closely related to the stochastic Taylor expansion, or series 
expansion of exponential maps on a noncommutative algebra. 

The actual simulation is carried out using (quasi) Monte Carlo methods. 



That is, one computes jfY^fLifi^T^'^ix)) where x!p^'^{x), i = 1,...,N, 
denotes N i.i.d. copies of x!^\x). Therefore, using the law of large numbers, 
the final error ^Ef=i f{X^^'\x)) - E[f{XT{x))] is of the order 0(^ + 
^). Then the optimal asymptotic choice of n is 0{n^) = 0{-\/N). 



The goal of the present article is twofold. First, we introduce a general 
framework to study weak approximation problems from the standpoint of 
operator (semigroup) expansions. That is given two processes that have 
equal semigroup expansions up to some order lead after composition to 
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two processes that are closed in law. This goal is not new. In fact, using 
PDE techniques, Milshtein and Talay between others proved various weak 
approximation results. Although our proof is essentially the same it gives a 
new viewpoint that will help in defining new approximation schemes. 

The next idea is to decompose the generator associated with (1.1) in (say) 
d + 2 components where each component is associated with each component 
of the driving process (the whole Levy process is considered as one compo- 
nent). Then we prove that if each of these components is approximated with 
an error of order m + 1 then the composition gives an error of order m. In the 
particular case that each component can be characterized as the semigroup 
of a flow-type process then the composition leads to a composition-type 
approximation scheme. 

Secondly, using the above strategy we provide approximations for solu- 
tions of (1.1). In particular, our approximations are valid for infinite activity 
Levy processes Y . We prove that in fact, if one uses the Asmussen-Rosihski 
idea of approximating the jumps of size smaller than e with a Brownian 
motion and we only simulate one jump of size bigger than e per each time 
interval in the approximation is enough to provide a first-order approxima- 
tion procedure. Furthermore we give the necessary estimate to determine e 
as a function of n. For this approximation, we found it better to decompose 
the generator in d -|- 4 components. 

This paper is organized as follows. In Section 2, we introduce the main 
example and the goal for the first part of this article in explicit mathematical 
terms. The general framework is introduced in Section 3. In Section 4 we give 
the results of convergence rates of numerical discretization schemes in the 
general framework. In Section 5, we give a general result that states how to 
recombine the approximations to coordinate processes in order to approxi- 
mate the semigroup associated to (1.1). Finally, in Section 6 we approximate 
each coordinate process and in particular, we define approximation schemes 
for Levy driven SDEs. 

2. Weak approximation problem. In order to better understand the ab- 
stract formulation in Section 3, we introduce here our main example. Let (1^) 
be a d-dimensional Levy process characterized by Levy-Khintchin formula 



(2.1) 



= e^pt(^^{e,b) - + J^^ie"^'^'^ - 1 - r(2/)))i/(d2/)) , 



where b G R'*, c G R'^ ® R'^ (symmetric, semi-positive definite) and z/ is a 
Borel measure on Rg := R'^ \ {0} satisfying that 



(2.2) 
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This measure v is called the Levy measure. It is well known that (2.2) 
implies that Yt G HpM-^^^ for all We also recall that r is a truncation 
function [e.g., r(y) =yl||j,|<i}, the constant b and r depend on each other]. 
The triplet (b, c, u) is called the Levy triplet. 

The Levy driven stochastic differential equation is given by 



Xt{x)=x+ [ %{XUx))ds+ f V{Xs-{x))dBs 

(2.3) 7 ^° 

+ / h{Xs-{x))dYs 
Jo 

with smooth coefficients Vq-.K^ ,V = {Vi, . . . ,Vd),h:R^ 0K'^ 

whose derivatives of any order (> 1) are bounded. Here Bt and Yt are inde- 
pendent d-dimensional standard Brownian motion and 1^ is a d-dimensional 
Levy process associated with the Levy triplet (6, 0,i/) satisfying condition 
(2.2). Using general semimartingale theory (see [17]) we have that the above 

equation has a unique solution. We define Vq :=Vq — ^ J2i=i J2jLi If"^^''^- 
Then (2.3) can be rewritten in the following Stratonovich form: 

Xt{x)=x + y2 fvi{Xs^{x))odBl+ f\{Xs^{x))dYs, 
■^Jo Jo 

where Bf = t. 

Before introducing the general framework of approximation, let us explain 
in mathematical terms the goal in this article. Our main example corre- 
sponds to the approximation of the semigroup Pt defined as the semigroup 
associated to the Markov process Xt'. 

Ptf{x) = E[f{Xt{x))], 

where / : ^ R is a continuous function with polynomial growth at in- 
finity. 

Let Qt = Qt be an operator such that the semigroup property is satisfied 
in {kT/n; k = 0, . . . ,n}. Assume that Qt approximates Pt in the sense that 
it satisfies the local error estimate {Pt — Qt)f{x) = 0{t"^~^^). Then using the 
semigroup property of both Pt and {QkT/n)^ notice that 

Prfix) - {QT/nTfix) = Y.(QT/n)HPT/n " Qr/n)^r-(fc+l)/nT/(^) • 

k=0 

Therefore if we have good norm estimates of {Qt/u)'^ and PT-(k+i)/nT in a 
sense to be defined later (in particular the norm estimates have to be inde- 
pendent of n) then we can expect that {Qt/u)^ is an approximation of order 
m to Pt- Finally in order to be able to perform Monte Carlo simulations we 
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assume that Q has a stochastic representation. That is, there exists a stochas- 
tic process M = Mt{x) starting at x such that Qtf{x) = £'[/(Mj(x))]. Then 
clearly, we have the following representation. 

Qrfix) = iQT/n)Vix) = i^[/(Mi/„ o . . . o M^/„(x))], 

where are independent copies of Mx/n and o is defined as (Mj o 

Mi)ix) ■.= MiiMlix)). 

The above ideas are well known and have been already used to achieve 
proofs of weak convergence (for historical references, see [9]). Nevertheless, it 
seems this is the first time it appears in this general framework. For example, 
if we take Mtix) := x + VQix)t + Vix)Bt + /i(x)yj for d = 1, one obtains the 
Euler-Maruyama scheme. 

Next to further simplify the procedure to obtain approximations we write 
the operator Pt as a composition of d + 2 operators as follows. First de- 
fine the following stochastic processes Xj^i(x), i = 0, . . . ,d + 1, usually called 
coordinate processes, which are the unique solutions of 

Xo^tix) = x+ f VoiXo^six))ds, 
Jo 

Xi^tix) = x+ f ViiXi^six)) o dB\, l<i<d, 
Jo 

Xd+i,tix) = x+ [ hiXd.+i,s-ix))dYs. 
Jo 

Then we define 

(2.4) Qi,tfix):=E[fiXi,tix))] 

for continuous function / : — > R with polynomial growth at infinity. 

For notational convenience we identify a smooth function V : R^ — > R'^ 
with a smooth vector field X^i^i^^'^^ R^. Let us define integro- 
differential operators Lj acting on by 

Lofix):=iVof)ix), Lifix):=\iV^f)ix), l<i<d, 

(2.5) Ld+ifix):=Vfix)hix)h 

+ j ifix + hix)y) - fix) - Vfix)hix)Tiy))uidy). 

It is well known that L := Yli=o generator of X and similarly Li 

is the generator of Xi^f Also e*^ := Pt and e''^' := Qi^t, respectively, where 
we consider these expressions as exponential maps on a noncommutative 
algebra. One notices that these operators have the form 

(2.6) e*^ = ^^L'= + 0(t'"+^), 

k=o 
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(2.7) e*^» = ^^Lf + 0(t 



To approximate e*^ , we would like to find some combination of operators 
satisfying 



k 



(2.8) e*^ - ijc^^'^^^'^ ■ ■ ■ = 0{t 



with some tjj > 0, G {Lq, Li, . . . , -L^+i} and weights {^j} C [0,1] with 
Sj=i'?i = 1- This will correspond to an jnth-order discretization scheme. 

To find such schemes, one can perform formal Taylor expansions for e*^ in 
each of the terms in (2.8). We remark that the result (2.8) will follow directly 
from (2.6) and (2.7) independent of the specific form of the decomposition 
L := X]i=o ^i- This algebraic calculation has lead to the introduction of the 
following approximation schemes: 

Ninomiya-Victoir (a): 

(2.9) ^e*^^^°e*"^^ • • • e*^''-+^e^^'^^° + le*^^"^°e*'^''+^ • • • ^^^^1"^^^ . 
Ninomiya-Victoir (b): 

(2.10) ^e^Lo^tLi__^tLa+i _^ igtid+i...gtiigtio_ 

Splitting method: 

(2.11) 6*/^-^° • • • e^/'^^'i^'tLa+i^t/2La . . . gV2io_ 

The semigroups generated by these operators have a probabilistic represen- 
tation. For example, Ninomiya-Victoir (a) corresponds to 

+ '^l/2<U^0,t/2 ° Xl,t ■ ■ ■ Xd+l,t O XQ^t/2{x), 

where ?7 is a uniform random variable taking values in [0, 1], independent of 
Xi^f However, since a closed-form solution Xi^t is not always available, one 
has to replace Xi t with other approximations of order m + 1 so that the final 
approximation result remains unchanged. Nevertheless the fact that there is 
only one driving process simplifies this task. This problem will be discussed 
in Section 5. 
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3. Preliminaries. 

3.1. Notation and assumptions. In this section, we consider a general 
framework for weak approximations following the arguments in Section 2, 
without using the specific form of the operator. We first define the following 
functional spaces. 

• C™ = C^(R^): the set of C"^ functions / : ^ R such that for each 
multi-index a with < |a| < m, \d"f{x)\ < C(a)(l + for some posi- 
tive constant C{a). 

We also let Cp = Cp. Let us define a norm on by 

ll/llc™ := inf{C7 > : \d^fix)\ < C(l + |x|^'), < |a| < m,x G R^}, 

where we denote \a\ := ai -|- • • • -|- on for a = (ai, . . . , qtv) G 

• Cj^''"([0, T] X R^): the set of functions / : [0, T] x R^ ^ R such that s ^ 
f{s,x) is continuous differentiable for all x G R^ and satisfies that f{s, •), 
dj{s,-) e with sup,g[o,T](ll/(5,-)lb;r + \\dsf{s,-)\\c^) < oo. 

From now on, we denote by Qt : Up>o ~^ Up>o ^pl-^^^) ^ linear 

operator for <t <T. 

Assumption (Mo)- If / g Cp with p>2, then Qtf G Cp and 

sup ||Qi/||c,</^||/||c, 
te[o,T] 

for some constant K > independent of n. Futhermore, we assume < 
Qtf{x) < Qtg{x) whenever 0< f <g. 

We now introduce two assumptions which are highly related to the con- 
vergence rate of approximation schemes. 

Assumption (M). Qt satisfies {Mq), and for each fp{x) := \x\'^p (p G 
N), 

(3.1) Qtfpix) < (1 + Kt)fp{x) + K't 

for some constant K = K{T,p), K' = K'{T,p) > 0. 

For m G N, 6m-[0,T] — > R+ denotes a decreasing function which satisfies 

1- ^m{t) „ 
lim sup r = U. 

Usually, we have 6m{t) = t^. 
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Assumption TZ{m,6m)- For each p>2, there exists a constant q = 
q{'m,p) > p and Hnear operators : Cp'' — > Cp-|_2fc (A; = 0, 1, . . . , m) such that 

(A) : For every / € Cp^"^ ^"^^ with 1 < m' < m, the operator satisfies 

m' 

(3.2) gj(x) = ^(efc/)(x)t'^ + (Err^'V)(x), [0,T], 

fc=0 

where ErrJ f £ Cq, and satisfies the following condition: 

(B) : If / G C;^" with m" > 2k, then e^f G C^'a^^^' and there exists a 
constant = K{T, m) > such that 

(3.3) ||efc/|| „//_2fe < /^||/|Um", A; = 0,l,...,m. 

Furthermore if / € C^" with m" > 2m' + 2, 

IIF (-')fll </^*'"''''ll^llcr"' if"^'<"^> 
" * \i^t5„,(t)||/||^^„", if m' = m 

for all < t < T. 

(C) : For every < /c < m and j > 2/c + 2, if / € CjJJ([0, T] x R^), then 

In order to compare the finite power expansions of different operators, we 
introduce the following notation: 



J<m{Qt) '■= X! ^^^k, 

k=0 

J<m{Qt) is a linear operator, which is related to the series expansion of 
1 1— > e*^' (cf. Proposition A. 6). The following lemma comprises some basic 
properties related to the above definition. The proof is straightforward. 

Lemma 3.1. The following properties are satisfied: 
7^(m + l,5^+l)^7^(m,^"^), 

n{m,5m) =^ TZim,6m), 
whenever 5m{t) < K6m{t) and limsup^_^Q_|_ 6m =0. 

(i) Let {£,i}i<i<e be deterministic positive constants with X^i'^j = 1) '^''^d 
assume (M. ) for Q^^ (i = 1, ...,£) . Then J2i=i CiQ^t^ o.^so satisfies (M.)- 

(ii) Let {^i}i<i<£ C R and assume TZ{m,5m) for q[^^ {i = 1, ...,£) . Then 
J2i=i^iQi^^ a/so satisfies TZ{m,6m)- 
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4. Weak rate of convergence. In this section, we prove the rate of con- 
vergence for the approximating operator Q under the Assumptions (A^), 
TZ{rn,5m)- Throughout this section, we assume the following assumption: 

Assumption {M p). For all / G then P./ G C^p^^"^ and furthermore 
the following two properties are satisfied for some positive constant C: 
1- SUPtg[o,T] ll^i/llc- <C||/||c-, 
2. \\{P^-P^)f\\c^<C\t-s\\\f\\c^. 

Theorem 4.1. Assume (M) andTZ{m,6m) for Pt andQt with J<m{Pt- 
Qt) = 0. Then for any f G Cp^™"*"^^ , there exists a constant K = K{T, x) > 
such that 

(4.1) |PT/(x)-(QT/n)"/(x)|<i^5™(^)||/||^2(,„+i). 

For the proof, we need the following lemma. 

Lemma 4.1. Under Assumption (M), the operators Pt and Qt satisfy 
sup rnax ((Pr/n)'' + iQT/n)'')f{x) < oo 

n 0<k<n 

for any positive function f £ Cp with p>0. 

Proof. Without loss of generality we do the proof for Q. Let fp{x) = 
|3.|2p £qj. p g By the Assumption (M), we have 

(QT/n)'7p(^) = {QT/n)''-\QT/nfp){x) 

<fl + -)(gT/n)'-Vp(x) + - 

with some constant C,C' independent of t,x,k,n. Since (1 + < e'^ , one 
proves by induction that 

sup max 

n 0<k<n 

This completes the proof. □ 

Proof of Theorem 4.1. Let / G Cp^™+^\ Using the semigroup prop- 
erty and Assumption TZ{m,5m), we have 

n-1 

Prfix) - (QT/nTfix) = Y.iQT/n)\PT/n " Qt /n)PT^{k+l)/nT f {x) 

k=0 
n-1 

= T.(QT/n)\^rr^^iPT-ik+l)/nTf){x), 
k=0 
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where Err|™^ is the error term of {P — Q) defined in (3.2). 
We obtain from Assumptions TZ{m,5m) and {Mp) 

|(Erri?;iPr-(fc+i)/nT/)(x)| < + m\PT^ik+i)/nTf\\c^i^+^) 

<^5™f-)(l + |xr)||/|L2(„.+i) 
n \n J 

and hence Lemma 4.1 leads to 

\{QT/n)H'^^4tPT-ik+l)/nTf){x)\ 

< (-) ii/iL2(™+i) (Qt/„)'(i + m 

< —Smi — ) ||/|L2(m + l) 

for some constant K = K{T, x). This completes the proof. □ 

The following theorem is an extension of Theorem 4.1, and is analogous 
to Talay and Tubaro [19], Theorem 1. 

Theorem 4.2. Assume (M ) and Tl{m + l,5m+i) for Qt with J<m{Pt - 
Qt) = 0. Then for each f G Cp^^~^'^\ we have 

(4.2) PTf{x) - {Qr/nTfix) = ^ + o[[^y^\ (^)) , 
where K = T^ PsJm+i{P - Q)PT-sf{x) ds. 

Proof. We start by noting that as in the proof of Theorem 4.1, 

{Px/n - QT/n)PT-sf{x) 

= j Jm+i{P - Q)PT-sf{x) + iErr^^^^'> PT-sf){x) 
and therefore, 

Prfix) - (QT/nTfix) 

= - T^iGr/rd^Jm+liP - Q)PT-ik+l)/nTf{x) 

^ ^ fc=0 

+ o(w(^)). 
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Now applying the proof of Theorem 4.1 (for m = 1) to Jm+i{P — Q) ^ 

PT-{k+l)/nTf G C^+2(m+l)' obtain 

l((QT/n)' - PkT/n)J.^+l{P-Q)PT-{k+l)/nTf{x)\ 

< ^^^^^\\Jm+l{P - Q)PT^(^k+l)/nTf\\cl^^^^^^^ 



< 



C2{T,x) 



n 



c 



2{m+3) • 



Next, we have by hypothesis (A4p), 

\PkT/nJm+l{P - Q)PT~(k+l)/nTf{x) 

- P{k+l)/nTJm+l{P - Q)PT-{k+l)/nTf{x)\ 

= \{I - PT/n)PkT/nJm+l{P - Q)PT^(k+l)/nTf{x) 

' ^ ^PkT/nJm+l{P — Q)PT-(k+l)/nTf\\c'-^ 



< 



< 



n 

Ci{T,x) 



p+2{m + l) 



-,2{m+3) • 



Using Lemmas A.l and A. 2 in the Appendix and Jm+i{P — Q)PT-sf{x) £ 

■<1,2 

■'p+2(m+2)' 



n-1 



A;=0 



~ X! ^{k+l)/nTJm+l{P - Q) Pt - [k+l) / nT f 

- r PsJm+l{P-Q)PT-sf{x)ds 

Jo 



< 



CiT,f,x) 



n 



As a result, taking K = T"' PsJm.+i{P-Q)PT-sf{x) ds, we conclude that 



Prfix) - (QT/nTfix) = 

This concludes the proof. □ 



K 



+ 



n" 



n 



m+l 



y 6, 



m+l 



5. Algebraic approximations of semigroup operators using coordinate op- 
erators. Throughout this section, we assume that Pt, t € [0,r], is a semi- 
group that satisfies {^A), (Aip) and lZ{m,5m)- Furthermore we suppose 
that 

m 



P 
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with ej = (X]f=o satisfying the properties stated in TZ{m, 5m)- Similarly, 
we assume that Qi^t '■ Up>o C'p(R^) — > Up>o C'p(R^), i = 0, . . . , d + 1, be a 
sequence of operators such that they satisfy (A^), (A^ p) and 7^(m, f^^) with 

t J' • 

J<rr^(Q^,^)=/ + E-^^ 

ni=i cii := (21(22 • • • Of denotes a noncommutative product. 



Theorem 5.1. Assume m = 2. That is, (M.) andTZ{2,52) are satisfied 
for Qit (i = 0,1, . . . ,d + l) . Then all the following operators satisfy (M ) and 
n{2,62): 

N-V (a); QI"^ = ^Qo,t/2 nf=l Qi,tQo,t/2 + ^Qo,t/2 U.i=l Qd+2-i,tQo,t/2- 

N-V (b); = 1 nf+o^ Q^,t + k Uto Qd+l-^,t. 

Splitting: Q^"''^ = Qo,t/2 • • • Qd,t/2Qd+i,tQd,t/2 • • • Qo,t/2- 

Moreover, we have J<2{Qt^) = J<2{Q?) = J<2{QI'^^) = ELom^*'- 
particular, the above schemes define a second-order approximation scheme. 



The proof of Theorem 5.1 is an application of Theorem 4.1. The condi- 
tions follow from the next lemma, together with an algebraic calculation as 
pointed out at the end of Section 2. 

This theorem can also be stated for third-order approximation schemes. 

Lemma 5.1. Let Q\ and : Up>o C'p(R-^) ^ Up>o C'p(I^^) two lin- 
ear operators and let QjQf be the composite operator. Then: 

(i) // (M ) holds for Q\, Ql, then it also holds for QjQ^. 

(ii) If TZ{m,5m) holds for Qj, Ql, then it also holds for Q\Q1. 

Proof, (i) is obvious. We now prove (ii). Let m' < m. We have by 
hypothesis that 

m' 

Qlfix) = T.iJkQlf){x)t'^ + (Errl'"' 
fc=0 

m' 

Qlf{x) = Y^{J,Qlf){x)t^ + (Err^ 

fc=0 

for / G Cp^^ ^'^\ P > 2. Furthermore there exists q = q{m,p) > such that 
Errl™''^V,Err^'''VeC,.Now we prove (A)— (C) in the definition of T^i^ui^ ^m) • 
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(A): Note that for / G Cp^™'+^^(R^), 

QlQ^fix) = Ql (j2{JkQV)i^)t' + (ErrJ"' . 

\k=0 I 

Since JkQtf ^ ^1+21^^^ ^ Qli^kQlf) can be written as 

m' —k 

{Q]{JkQimx)= iJ,QHJkQimx)t' + {ETTt'-'^'\hQif)ix). 

e=o 

As a result, we have 

m' m' — k 

(m',1,2) 



where 



(5.1) (ErrJ-' = (Q.^ErrJ™' + Y.iM"''~'''^JkQ'tf)ix)t\ 

k=0 

im' 1 2) 

We obtain from the properties of the error terms that ErrJ ' / E Cqi for 
some q' = q'{m,p) > q. 

(B): For /EC™ with m" > 2(m' + 1), we can derive for k + i<m' 

\\JeQt {'JkQtf )\\r<"^"-2(k+t) < -fCill JfcQj /|Lm"-2fc < -fsr2||/|lc<m" 
and by (5.1), 

llErrf'^'^VlIc , < i^sllErrl-' + K4E4""''hoQjf\\c , 



+ K5j2\\JkQ'if\\c'-:'-2.t"''^' 



k=i 

ri^t'"'+i||/|L„», if m'<m, 
< % 

\^Kt6m{t)\\f\\Qn," , iim' = m. 
Finally, the proof of (C) is straightforward. □ 

Proof of Theorem 5.1. Using this lemma, we end the proof, calcu- 
lating J<m for each numerical discretization scheme. For instance, in the 
case of N-V (b) [i.e., (2.10)], we obtain 

'^<2 9 n + 9 n Qd+l-i,t 

i=0 ^ i=0 J 
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1 C^^^ \ 1 C^^^ \ 

-J<2 \J[ J<2iQit)j + ^J<2 \\[ J<2{Qd+l-i,t)j 

1 / d+l .2 d+l \ 

\ j=l j=l i<j / 

-I / d+l ,2 \ 

\ i=l i=l i>i / 



j<2{Pt: 



□ 



Another idea to construct higher-order schemes is to use local Romberg 
extrapolation. In order to do this we need to weaken the assumption {^j} C 
[0, 1]. This is done in the next theorem. 

Theorem 5.2. Let m = l or 2. Assume (M) and TZ{2m,t^^) for Pt 
and qI*' (i = 1, . . . , and (Aip) for Pt . Furthermore, we assume: 

(1) J<2m{Pt - YA=iiiQt) = for some real numbers {^i}j=i,.,./ with 

(2) There exists a constant q = q{m,p) > such that for every f € with 



m'>2(m + l), (P,-Ql*l)/GCg 



m'-2(m+l) 



anc 



sup \\{Pt-Qr)f\\^,..'^2i^+i,<CT\\f\\c^'r 



m+l 



te[o,T] 



Then we have for any f G Cp^'^'^^^ , 



PTf{x)-Y.UQ%nTfi^) 



i=l 



< 



CiTJ,x) 



n 



2m 



Proof. We first remark that the operator X]i=i dQt longer satisfies 

the semigroup property, that is, X^Li CiiQx/J'^ ^ iJ2i=i iiO^T/nT- ^^^^ 
proof is nontrivial. 
Note that for / G 

8 ■.= PTf{x)-j^UQ^^ rf{x)=j2UPT - (gP/J")/(x). 



i=l 



i=l 
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Using the semigroup property of Pt and Qi./^j', we have 



k/nT' 

I n-1 

j=l fc=0 

i n-1 

= PkT/n{PT/n " QT/n)PT-{k+l)/nTf{x) 

i=l k=0 

I n-1 

+ E^* E((QT/n)' - PkT/n){PT/n " Q^^ , J Pt ~ {k+l) / nT f {x) 
i=l k=0 

We expand (Qx/n)'' ~ PkT/n again, to obtain 

n-1 / ^ \ 

£=Y. (Pt/u)" Pt/u - E ^iQr/n PT-(k+l)/nTfix) 
k=0 \ i=l ) 

I n-lk-1 
+ Ee.EE(QP/n/(QP/n-^'T/n) 



i=l k=0 1=0 



X PT-{l+l)/nT{PT/n " Q^T / n) - {k+l) / nT f {x) . 

By the assumption (1), we have 

C7i(r,/,x) 



n-1 f ^ \ 

Y.{PT/nf Pr/n -Y.^^Q^T|n PT-(k+l)/nT f 
k=0 \ i=l ) 

Thus we end the proof by showing that 

I n-lfc-l 



< 



„2m 



i=\ k=0 1=0 



-{1+1)1 nT 



X 



{Pt/u - QT/n)PT-{k+l)/nTf{x) 



— j^2m 



Using here the assumption (2), we obtain 

IK^T/n ~ PT/n)PT-{l+l)/nT{PT/n " Q^T /J^T-ik+l) /nT fWc ^, 
- Z^^W^^T/n - QT/n)PT-{k+l)/nTf\\(j^(.m+l) 



< 



C'{T) 

2(m+l) II'' IICp 



n 



-^4{m + l) 
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and therefore 

I n-lfc-l 

-(l+l)/nT 

2=1 fe=0 /=0 



{Pr/n - QT/n)PT-{k+l)/nTf{x] 



< g g C2{TJ,X) ^ C2{T, f, X) 



„2(m+l) — ^2m 
fe=0 1=0 ^ ^ 

This completes the proof. □ 

Example 5.2. It is known that the Ninomiya-Victoir scheme 

/t d+1 1 d+1 \" 

£er/(2n)Lo -Q gT/nL,gT/(2n)Lo _^ igT/(2n)Lo "Q /nLa+2-^ J / {2n)Lo 

i=i ^ i=i / 

is of order 2 [m = 2, 52{t) = in Theorem 4.1]. By Theorem 5.2, the fohowing 
modified Ninomiya-Victoir scheme 

1 / d+1 \ " 1 / 'i+i \ 

i j^gT/(2n)Lo J-J gT/nL,gT/(2n)Lo j _^ £ / T/(2n)Lo "Q J /nLi+^-i J / {2n)Lo ^ 

is also of order 2. 

Example 5.3. Fujiwara [6] gives a proof of a similar version of the 
above theorem and some examples of fourth and sixth order. We introduce 
the examples of fourth order: 

In order to complete the approximation procedure through (quasi) Monte 
Carlo methods we need to find a stochastic characterization of the operators 

Qi,t- 

Definition 5.4. Given a stochastic process Yt{x) G Pl^^^L^, we say 
that Y is the stochastic characterization of the linear operator Qt if Qtf{x) = 
E[f(Yt{x))\ for / G Up>o C'p- In such as case we use the notation Qt = Qj ■ 

Remark 5.5. Given the operators Qf^ = 1, . . . ,£) and the determin- 
istic positive weights {6}i<i<^ with Yl\=i = 1- Let ?7 be a uniform random 
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variable on [0, 1] independent of and define Z := X^Li l(S}=i 6' ^ ^ < 
E;=lCi)^^ Tlien 

gf/(x)^i?[/(Zt(x))]=^6Qf7(x). 

i=l 

Tlierefore by Lemma 3.1 if satisfy (Al) and Tl{m,6m) so does Qf . This 
property will be used repeatedly in what follows. 

6. Applications. From this section on, we discuss the application of the 
previous approximation results to the case of solutions of the sde (1.1). Prom 
the results in the Appendix (see Corollary A. 7), it is clear that the semi- 
group Ptf{x) := E[f {Xt{x))\ satisfies the hypotheses {M) and TZ{m,5m)- 
We define various approximations generated via a stochastic process Xi 
with corresponding operator Qf^ (i = 0, 1, . . . , d + 1). 

Due to the previous results and in particular, Theorem 5.1, we see that is 
enough to verify local conditions on the approximation operators to conclude 
global properties of approximation. In particular, we only need to verify 
that the operator associated with Xi (the approximation to the coordinate 

process) satisfies (A^) and TZ{m,5m) and J<m{Qt^') = I + J2T=i7\Li for 
some m > 2 for Li given by (2.5). This is the goal in most of the applications 
in this section. 

Recall that the stochastic differential equation to be approximated is 

Xt{x)=x + y] fvi{Xs-{x))odBi+ f\{Xs-{x))dYs. 

In each of the following sections we consider different approximation pro- 
cesses for the coordinate processes Xi^f In each section, the notation for 
the approximating process is always Xi^t- We hope that this does not raise 
confusion as the framework in each section is clear. 

6.1. Continuous diffusion component. 

(a) Explicit solution. Let V : — > be a smooth function satisfying 
the linear growth condition |V^(2;)| < C(l + |x|). The exponential map is 
defined as exp(y)x = zi{x) where z denotes the solution of the ordinary 
differential equation 

(6.1) ^^ = V{ztix)), zo{x) = x. 

The solution of the coordinate sde is obtained in the following proposition. 
The proof follows from Ito's formula. 
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Proposition 6.1. For i = 0,1, . . . ,d, the stochastic differential equation 

(6.2) Xi,t{x) = x+ f V,{XUx)) o dBi 

Jo 

has a unique solution given by 

Xi,t{x)=exp{BlVi)x. 

Xi^t{x) is called the ith. coordinate process and its semigroup is denoted 
by Q\. This is a trivial example of the approximation of e*^' , i = 0, 1, . . . , d, 
satisfying {M) and lZ{m,t'^). However, sometimes it is not easy to obtain 
the closed-form solution to the ODE (6.1). In those cases, we shall approx- 
imate exp(tV)x. Here we will do this with the Taylor expansion first and 
then the Runge-Kutta methods denoted by hm and Cm, respectively. 

(b) Taylor expansion. We first prove the following lemmas which help 
us to find the rate of convergence of the scheme to be defined later. The 
following lemma follows easily from Gronwall's lemma. 

Lemma 6.2. Let V he a smooth function which satisfies the linear growth 
condition. Then \ exp(ty)j;| < Ce^l*l(l + |x|) /or t G R, x G R^. 

From now on we denote by Cj : R^ — > R, the coordinate function Cj [x) = 
Xj for J = 1, . . . , A^. Furthermore, we also denote by V the vector field oper- 
ator defined from V. 



(6.3) /(exp(tV-)x) = ^ TiVl'fix) + / ^ ^y-+i/(exp(nV,)x) du 

t^n Jo ml 



Lemma 6.3. Let f G C^~^^. Then we have for i = 0,1, . . . ,d 

k=0 

and 

f ^^^^Vl^+^f{eMnVi)x) du < C„||/||^™+ie^l*l (1 + \x\P+^+')\tr+^ 
Jq ml 

for all t G R. 

Proof. Assertion (6.3) follows application of Taylor expansion to the 
function f{exp{tV)x) around t = 0. Next, as \ Vl^^'^f{x)\ < C{1 + \x\p+"'+^) , 
we obtain from Lemma 6.2, 



* {t - uY 



V^+^f{exp{uV)x)du 



rnl 

\t\ 



<Cm\\f\\r^+i I |trCe^l"l(l + |xr"^+^)du 
p Jo 

<c'j\f\i.^+.e^\\i+\xr^+^)\tr+\ 
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□ 

Based on this lemma, we define the approximation to the solution of the 
coordinate equation (6.2) as follows 

64 (t, V)x = ^ j = l,...,N. 

k=o ^■ 

Define 

Xi^tix) = b2m+i{Bl, Vi)x for z = 0, . . . , 
Then we have the following approximation result. 

Proposition 6.4. (i) For every p> I, 

\\Xi,t{x) - X,,t{x)\\LP < C{p,m,T){l + 
(ii) Let f £ Cp. Then we have 

E[\f{X.A^)) - f{X.,ixm<C{7nX)\\f\\c-^^ 

Proof, (i) Apply Proposition 6.1 and Lemma 6.3 with / = Cj. Then we 
have 

||X,,i(x)-X,,i(x)|Up<i?[|C^e^^l^'l(l + |x|2{™+i))|i?,|2(-+i)ni/^ 
<C(l + |2;|2('"+i))t"^+i 

for some constant C = C{p,m,T). 

(ii) We first apply the mean value theorem to obtain 

E[\f{X,,t{x))-f{X.,t{m 

< + \eXi,t{x) + (1 - 9)x,4x)\p\\L4Xi,t{^) - Xi,t{x)\\L^ 

< c\\f\\ci\\i + \x,,tixW + \x,,t{x)nL4i + |x|2{-+i))r+^ 

We see by Lemma 6.2 that 

sup 111 + \Xi,t{x)\P + \Xi,t{x)\P\\L2 < C'{1 + \x\P) 
tG[0,T] 

from which the proof follows. □ 

As a result of this proposition we can see that TZ{m,t'"^) holds for the 
operators associated with bmit,Vo)x and 62m+i(-B^, Vi)x, l<i<d. Indeed, 
we have for m' <m, 

E[f{X,,t{x))] = E[f{X,,t{x))] + E[fiX,,tix)) - f{X,,t{x))] 

= Ery^'/(^) + (^r7)(^), 

fc=o 
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where 

{Er'f){x) := {Errt'^f){x)+E[fiXuix)) - /(^.^(x))] 

and (Err^"^ is defined through the residue appearing in Proposition 
A. 6, using Li and Qi instead of L and P. Furthermore, using Proposi- 
tion 6.4(ii), we have that the error term E"™ satisfies (B) in Assumption 
7^(m,^'"). 

It remains to prove that (M.) holds for ^(x). For the proof, we need an 
additional growth condition for the vector field V^. 



Proposition 6.5. Assume that {V^Cj) {2<k<m, <i <d, I < j < 
N) satisfies the linear growth condition then (M) holds for Xi^t{x), i = 
0,...,d. 



Proof. The Assumption (A^o) follows from the smoothness and the 
linear growth property of V^ej. We only prove the moment condition (3.1) 
for Xi^tix) i = 1, . . . ,d. Consider the multiplication {p E N) 



k=0 



2p 



+ BlV,{x) + Y^^-^{Vte,){x) 



k=2 



2p 



Taking into account that £^[(5^)^^+^] =0, A; G N. Then by the assumption, 
we obtain the result. □ 



Therefore we obtain the main result. 



Theorem 6.1. Assume that (V/'ej) {2<k<m, 0<i<d, I < j < N) 
satisfies the linear growth condition. Let Xi^t{x) be defined by 

X,4x) = b2m+i{Bl,Vi)x= E TrAV^I){x) / lodBi^o...odBl. 

Denote by Qf * the semigroup associated with Xi^t{x). Then ' satisfies 
(M ) and n{m, f^). Furthermore J<Mf') = 1 + TJLi %L{. 

(c) Runge-Kutta methods. We say here that Cm is an s-stage explicit 
Runge-Kutta method of order m for the ODE (6.1) if it can be written in 
the form 

s 

(6.4) Cmit,V)x = x + tJ2PMt,y)x, 

1=1 
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where ki{t,V)x defined inductively by 
ki{t,V)x = V{x), 

ki{t,V)x = V ^x + tY^^aijkj{t,V)x^ , 2<i<s, 
and satisfies 

I exp(ty)x - c^(t, V)x\ < C^e^l*l(l + \xr+^)\tr+^ 

for some constants ((/?«, ajj)i<j<j<s). Runge-Kutta formulas of order less 
than or equal to 7 are well known. For details, see, for example. Butcher [4]. 

The following proposition can be shown by the same argument as in the 
proof of Proposition 6.4. 

Proposition 6.6 (Stochastic Runge-Kutta). (i) For every p> 1, 
(6.5) \\Xu{x) - C2m+i{BlVi)x\\LP < C{p,m,T)il + |x|2(-+i))t™+^ 
(ii) Let f Cp. Then we have 

E[\fiX,^tix))-f{c2.m+iiBiVi)x)\] 

(6.6) 

< C(m,r)||/|bi(l + |x|2('^+i))t'"+i. 

Next we show that (M) still holds for the Runge-Kutta schemes. 

Proposition 6.7. (M ) holds for Cm.{Bl, Vi)x, i = 0,...,d. 

Proof. We first note that for every 1 < j < s, there exists a function of 
the form pj = J2k=Q^jk\A'' such that 

\kj{t,V)x\ <pj{t){l + \x\). 

Assumption (A^o) follows from the smoothness and the linear growth prop- 
erty of Vi. We now prove (3.1). In the case i = 0, this is obvious by definition 
and the inequality (6.1). In the case 1 <i <d, observe that 

Crr.(t,V)x = x^ty^^iV{x)^ty^^i \ —V\x^Qty^ai^jk.j{t,V)x\dQ 

1=1 1=2 V 1=1 J 



=:x + tY^I5iV{x)+ D^{t,V)x. 
1=1 

Expanding multiplications and taking expectations, as in Proposition 6.5, 
we can show that the terms containing odd powers of Bl have expectation 
0. Finally, we obtain from the boundedness of dVi that 

\D,n{BiV,)x\<p{Bl){l + \x\), 
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where p = p{t) is of the form X]|=2 '^^fcl^l'^- Using this, we conclude the proof. 
□ 

Consequently, as in the Taylor scheme, TZ{m,t'^) and {M) hold for the 
operators associated with Cm(t, Vo)x and C2m+i{Bl,Vi)x, 1 <i <d. For more 
on this method, we refer the reader to [12]. 

(d) Minor extension. In the previous approximation, the assumption 
that Bt ~ A^(0, Id) can be weakened. In fact, we can use \/tZ instead of 
Bt where {Z'^)f^i are independent and 

P{Z' = ±V3) = l, P{Z' = 0) = l 
for each i = 1, . . . ,d. 

Proposition 6.8. Let Bt be a one-dimensional Brownian motion and 
Z be a R-valued random variable such that for all <k < m, 

E[{Z)'] = E[{B^)'] 

and 

E[exp{c\Z\)] < oo 
for any c> 0. Then, for every f G C^~^^ , 

|i?[/(exp(Sty)x)] - i?[/(c^(^/tZ,y)x)]| < C(m,r)(l + |xr™+i)t("^+i)/2_ 

6.2. Compound Poisson case. Suppose that is a compound Poisson 
process. That is, 

Nt 

1=1 

where (Nt) is a Poisson process with intensity A and (J,,) are i.i.d. R'^-valued 
random variables independent of (Nt) with Jj G np>i-^^- 
In this case 1^ is a Levy process with generator of the form 



{f{x + y)-f{x))v{dy), 



where r = 0, 6 = 0, z^(Ro) = A < oo and v{dy) = XP{Ji G dy). 
Then in this case 

(6.7) Xf+\x)=x+ f\{Xf+\x))dYs, te[0,T], 

Jo 
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which can be solved exphcitly. Indeed, let {Gi{x)) be defined by recursively 

Go = X, 

Gi = Gi-i + h{Gi-i)Ji. 

Then the solution can be written as Xf^^{x) = Gi^^{x). Define for fixed 
M G N, the approximation process X^+i^t = Gn^;^m{x)- This approximation 
requires the simulation of at most M jumps. In fact, the rate of convergence 
is fast as the following result shows (see also Mordecki et al. [14]). 

Theorem 6.2. Let M G N. Then the process GNtAAiix) satisfies (M) 
and Tl{M,t^^~'^) for arbitrary small k> 0. Furthermore J<MiQf'''^^) = I + 
^i=i j\^d+i- 

Proof. Note that for f e Cp 
Qf'^'f{x)-Q'}+'f{x) = E[f{GN,AM{x))]-E[f{GNAx))] 

where Tm := inf{t > : A'^t = M}. By the Holder inequality, 

|gf^+7(^)-Q^V(x)| 

< 2E 



2E 



sup |/(G^,(x))r/(^-i) 

0<t<T 

sup |/(G^,(x))r/(^"i) 



{7-l)/7 

P{Tm+i < tf'^ 



1 {7-l)/7 / /•< (As)^ - N 



LO<t<T 

<C(7,T)||/||c,(l + |xn(tA-i)(^-^+i)/^. 







Ml 



-Xe-^' ds 



K : = 



Take sufficiently small 7 > 1, then TZ{M,t^'^~'^) holds for Qf''^^ where 
(1 - 1/7)(M + 1) > 0. Finally, we show [M). Let fp{x) = (p g N) and 
7 < M . Then using the above calculation and Corollary A. 7, we have 

< (1 + Krt)fp{x) + K2t + IQf'+'Uix) - Qt'fp{x)\ 

< (1 + K^t)U{x) + Kit. □ 

6.3. Infinite activity case. In this section, we consider the SDE (2.3) 
under the conditions z^(Ro) = 00. Without loss of generality, we assume 
that c = 0. 
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(a) Ignoring small jumps. Define for e > the finite activity (i.e., drift 
+ compound Poisson) Levy process (Yf) with Levy triplet {b,0,u^) where 
the Levy measure is defined by 

(6.8) u'{E):=iy{En{y:\y\>£}), E e BiR!^). 

Consider the approximate coordinate SDE 

Xd+i,t{x) = x+ I h{Xd+i,s-{x))dY,^, 



result shows (see also Mordecki et al. [14]) whose generator is 



Ll'lJix) = Vf{x)h{x)b + I {fix + h{x)y) - fix) - V fixMxMy)^ {dy). 



Now we derive the error estimate for X^+i.t- 

Theorem 6.3. Assume that < e = e(t) <1 is chosen as to satisfy that 
a'^ie) := /|j^|<£ \y\'^i'idy) < Then we have that Qf''^^ satisfies (M ) and 

7^(M, t^'). Furthermore J<A/(Qf '+') = / + ET=i jA^v 

Proof. First, we remark that condition (A^o) follows from Proposition 
5.2 in [7]. We start by noting that from Proposition A. 6, we have 



QrV(x) - Q^+7(^ 



(6.9) 



fe=i 



+ 



'''^''iQi^\Ld^ir+'-Q^'-\L]ii,r^')fix)du. 



Ml 



Therefore the proof is achieved if we prove that 

|(L,+i - Ll'l,)fix)\ < C\\f\\c.il + \xr')t''+\ 

For the proof, we change here the representation of the Levy triplets of Yj 
and Yf as follows: 

ib,0,u),T^ ibe,0,lj),Te, 

(6,0,i/^),r^(6„0,i/^),r„ 
where Tsiy) =yl{|y|<e}- Then 



l(L,+i-4:f,)/(x)| 



(6.10) 



< 



Vfix)hix)iy - Teiy))ividy) - v'idy)) 



d 



1 - ^)^/(^ + ^hix)y) ddividy) - v'idy)) 



AN OPERATOR APPROACH FOR MC WEAK APPROXIMATIONS 



25 



We first obtain that for e > 0, 

{y-TM)Hdy)-u'{dy)) = Q 



since the support of the measure {v — u^) is {|y| < e}. Now we consider the 
second term of (6.10). We can immediately show that due to the polynomial 
growth property for /, 

"1 <P 



^^^f{x + eh{x)y) de{u{dy) - v'{dy)) 



<C\\f\\Ml + \xr^)a\e) 



and hence as a'^{e) < ^*^+^ one obtains that J<M{Qf^^^) = I + T,f=i 

and that satisfies {M) and n{M,t^^) follows as in the proof of Propo- 

sition 6.2. □ 

Using Theorem 5.1, we can incorporate the above approximating process 
^d+i,t to the whole approximation method. This will require to first simu- 
late the jump times of the approximating Levy process and then solving 
ODEs between these times. If the task is time consuming one can also sepa- 
rate the jump component from the drift component as indicated by Theorem 
5.1 (see also Section 6.4). The right size of e is determined by the condition 

(b) Approximation of small jumps. We apply here the Asmussen-Rosihski's 
approximation for small jumps of Levy processes. The idea is that the small 
jumps ignored in (6.8) are close to a Brownian motion with small variance 
(7^(e) (see details in [2]). 

Consider the new approximate SDE 

(6.11) Xd+i,t{x) = x+ [\{Xd+iAx))^\'^dWs+ f h{Xd+i,s-{x))dY^, 

Jo Jo 

where Wt is a new d-dimensional Brownian motion independent of Bt and 
Yf , and is the symmetric and semi-positive definite dx d matrix defined 
as 



(6.12) S,= / yy*i^{dy). 

J\y\<e 

We remark that is of the form AAA* , where A is an orthogonal matrix 
and A is the diagonal matrix with entries Ai, . . . , > (eigenvalues). Thus 

we use the notation = AA^/^. Since the above SDE is also driven by 
a jump-diffusion process, we can also simulate it using the second-order 
discretization schemes in Theorem 5.1. 

Now we prove that rate of convergence in this case is faster than in the 
case that we ignore completely the small jumps (see Theorem 6.3). 
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Theorem 6.4. Assume that < e = e{t) < 1 is chosen as to satisfy 
that /|y|<e lypi^(rfy) <t^^^^. Then we have that Qf"^^^ satisfies (M) and 

7^(M,^*0• Furthermore J<M{Qf'^') = I + Y.T=ijA^i- 

Proof. As before, condition (A^o) follows from Proposition 5.2 in [7]. 
The SDE corresponds to the generator 

L%J{x):=Vf{x)h{x)b+\Y.dk,if{x){h{x)T.,h*{x))k,i 

^ k,i 

+ I {fix + h{x)y) - fix) - Vfix)hix)riy)Xidy). 
Using this representation, we have for f ^C^, 

iLa^,-L%,)fix)= f f \l-e)^fix + 6hix)y)deiuidy)-u%dy)) 



^J29k,ifix)ihix)J:eh*ix)h 

k,l 

1 (1 - 6')2 d3 



2 de^ ^^"^ + ^^^""^^^ dOii^idy) - u'idy)). 
Hence we finish the proof as in the proof of Theorems 6.2 and 6.3. □ 



If we put all the pieces together, we have the following final result. Here 

t 



denote i = l,...,d,j = l,...,2n denote 2nd independent standard Brow- 



nian motions and B^-' = t. 

Theorem 6.5. Assume that Vq, V and h are infinitely differentiable 
functions with bounded derivatives with J-^dil A\y\P)uidy) < oo for allp G N. 

Define e = eiT,n) so that J^y^^^\y\'^iyidy) < i^)^ . Let Xl^ix) = c^iBl\Vi)x, 

1 = 0, . . . ,d, j = 1, . . . , 2n, 2n copies of the Runge-Kutta method of order 

2 as defined in (6.4) and j = l,...,2n, independent copies of 

the approximation defined in (6.11). Then the following schemes, X^"^ = 
o Y^~^ o • • • o Y^ix), are second-order discretization schemes: 
N~ V (a) .• Yi, ix) = f/, o o . . . o o X^^^^^^^ (x) + ( 1 - 

^i)^o,T/(2n) ° ^d+i,T/n ° " " " ° ^,T/n ° ^o,T/(2n)(^) ""^^re Uj is a Bernoulli 
r.v. with PiUj = 1) = 1/2, independent of everything else. 

N-V ih): Ylix) = V,X\^^^^i^ 0...0 Xl^i^ix) + (1 - U,) o • • • o 

^d+iT/n('^) ^^^'"c a Bernoulli r.v. with P{Uj = 1) = 1/2, independent 

of everything else. 
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Splitting: Yl{x) = ^^,r/(2n) ° ' " ' ° ^,T/{2n) ° ^d+i,T/n ° ^d,T/(2n) ° • • • ° 

^0,r/(2n)(^)- 

One can also write a similar result for higher-order schemes using Theorem 
5.2. 

6.4. Limiting the number of jumps per interval for approximations of in- 
finite activity Levy driven SDEs. In the previous two approximations al- 
though £ G (0, 1) may be relatively large compared with the interval size 
T/n, one still faces the possibility of having many jumps in the interval 
[0,T]. Therefore we introduce the idea used in Section 6.2. That is, we pro- 
pose another approximation that restricts the numbers of possible jumps to 
at most re. Throughout this section we assume that f\y\^i \y\i^{dy) < oo and 
without loss of generality, we assume that r(?/) = yl^y^^i. 

Then we decompose the operator 

L\^J{x):=Vf{x)h{x)(b- [ T{y)uidy)), 

V Je<\y\<l J 

Ll+if{x) := / (/(x + h{x)y) - f{x) - V f{x)h{x)T{y))u{dy), 
J\y\<e 

Ll+if{x) := I f{x + h{x)y) - f{x)u{dy). 
Je<\y\ 

The operator -L^^j^ can be easily approximated using any Runge-Kutta 
method for the ordinary differential equation 

Xj+i,, = x+(h- f r{y)u{dy)) f h{Xl^^^,) ds. 
V ''e<\y\<l / JO 

We denote by X^_^^ ^, the Euler scheme associated with this ordinary differ- 
ential equation. Therefore we only need to approximate L'^j^^ and L'^_^]^. 

Let I : R'^ be a localization function that may be used for im- 

portance sampling of the Levy measure. Let F^{dy) = Xj^l{y)l\y\<:^i'{dy) 

with A, = J\y\<J{y)i^idy). Let Y, ~ F,. Define xf''{x) = x + h{x)WtVX;, 
where W is a d-dimensional Wiener process with covariance matrix given 
by Sjj = l{Y^)~^YfYj which is independent of everything else. 

First we prove that X^'^{x) satisfies Assumption (Ai). 

Lemma 6.9. Assume that for p>2, sup^g^Q ^ x 
i^{dy) < oo, then Assumption (M ) is satisfied with 

E[\X%,{xW]<{l + Kt)\x\^ + K't. 
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Proof. Let f{x) = \x\^, p > 2. Using Ito's formula for p ^ 3 and an 
approximative argument in the case p = 3 (as in the proof of the Meyer-Ito 
formula) one obtains that 

(6.13) E[fiX^^'{x))]-f{x) 

(6.14) =^\,e' 



^i^T'l i^l-l^\X^'%x)r\h{x)Y^,X^'%x)f 

+ \X'^''{x)\P-^\h{x)Y'\^ds 



We use the Lipschitz property of h to obtain that 

\X^,'%x)\ = \x + hix)W,^,\ 

<{l + C\Ws\y/X'e){l + \x\). 

Then, we have 

\E[f{X^'%x))]-f{x)\ 

< Cpt{l + \xn f \y\\l + i\y\%y)-^Xet)^P-''y'Hdy). 

Lemma 6.10. Assume that for p> 2, 

Mp= sup / \y\^liyrHl + i\y\^liyr'X,t)^P-^^/^Mdy)<oo 

ee{0,l]-'\y\<e 

and J\y\<s \y\^i^{dy) < Ct then 

\E[f{X?^')] - fix) - tLl^J{x)\ < C{p)\\f\\c.{l + \xr')t\ 
That is, x'^'^{x) satisfies Assumption TZ{2,t^). 

Proof. Let f then applying Ito's formula, one gets 



□ 



E[f{X^^^)]=f{x) + ^E 



J2 d,jf{Xl'')hikhu{x)l{Y')-'Y^Yfds 



i,j,k,l 



fix) + 7; f Y] dijf{x)hikhii{x)ykyiu{dy) + Reix), 

2 J\v\<e 



where by Lemma 6.9, we have 

1^.(^)1 <^^ii/ib(i+i^r')*' 



X / \y\%y)-\l + i\y\^l{y)-'Xet)^P-'^/^)i^idy). 
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Furthermore 



Ldlif{x)-l; I V dijf{x)hikhii{x)ykyiv{dy) 

y2 [ [ {dijf{x + ah{x)y)-dijf{x))adahikhii{x)ykyiv{dy). 

~f.J\v\<sJQ 



Therefore 



- \ I dijf{x)hikhii{x)ykyiu{dy) 

<C\\f\\cfM + \xr') [ \y\Mdy)- 
J\y\<e 

This finishes the proof. □ 

In the particular case that l{y) =y^ , r = 2, the above scheme corresponds 
to a Asmussen-Rosiiiski type approach. 

The approximation for L^_^_i is defined as follows. Let 

Ge,i{dy) = C^il{y)l\y\y^u{dy), 
Ce,i = I l{y)v{dy) 



and let Z^'' ~ G^^; and let S'^'^ be a Bernoulli random variable independent 
of Z^'K Then consider the following two cases. If 5^'' = define Xf\x) = 
X, otherwise Xf'^{x) = x -\- h{x)l{Z^'^)~^ Z^'K Then we have the following 
results. 

Lemma 6.11. Assume that for p> 2, sxiY>i,^(^Qi^J]^y]^-^^l{y)~P\yY''^^v{dy) < 
oo and C~Ip[S^ = 1] < Ct then Assumption (Ai ) is satisfied with 

E[\x'4,ix)\P]<{l + Kt)\x\P + K't. 

Proof. The result follows clearly from {f{x) = 
P[S' = l]\E[f{x + h{x)l{Z^''r'Z^'^) - f{x)]\ 

= CrlP[S' = 1] [ {fix + h{x)liy)-'y) - f{x))l{y)u[dy) 



< Ct{l + IxD (l+ I l{yr^\y\^+^v{dy)] 
V J\y\>e J 



□ 
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Lemma 6.12. Assume that for f £Cp, we have that /|j^|>£ ^ 
1) + \yK^\l{y)~^ - iK^i^idy) < C and \C^^P[S^'^ = l]-t\< Ct^ the 

\E[f{Xf')] - fix) - tLl^J{x)\ < Ct'WfWciil + \xr'). 



len 



Proof. As before let f ^ Cp then 



E[f{xf'')] = fix) + E[f{x + h{x)l{Z^'r'Z'') - /(^); 

= fix) + / (fix + h{x)liy)-'y) - f{x))l[y)u{dy) 
■J\y\>e 

xC~/P[5^'' = l]. 

Then we clearly have that 

\E[f{xf')]-f{x)-tLl^J{x)\ 

fi 



<t 



\y\>eJid 



Y^{dif{x + ah{x)l{y)~^y) 



dif{x + ah{x)y)) da h{x)yv{dy) 



^C;lP[S''' = l]+ I f{x + h{x)y)-f{x)v{dy) 

x\C;}P[S'^' = l]-t\ 

<C\\f\\ci{l + \xr')t\ 
This finishes the proof. □ 



Using the previous results we can propose various schemes of approxima- 
tion of order 1 as in Theorem 6.5. We state the simplest type of approxima- 
tion. 

Theorem 6.6. Assume that Vq, V and h are infinitely differentiable 
functions with bounded derivatives with J-^d{l A\y\^)i'{dy) < oo for allp G N. 

Define e = £{T,n) so that the conditions on Lemmas 6.9, 6.10, 6.11 and 
6.12 are satisfied for t = T/n and for appropriate localization functions. Let 
Xl^{x),i = 0, . . . , j = 1, . . . ,n, n copies of the Euler-Maruyama method 
for Xi^tix). 

Also, let ^l^^ix/n' * ~ 1' 2' 3? i ~ 1' ■ • ■ ' ^? be n independent copies of the 
schemes defined above. Then the following scheme, x!^'^ = o Y^^"^ o • • • o 

Y^{x), Yl{x) = X^^i^ o . . . o o Xl'll^^^ ° ^I'+It/u ° ^d+l,T/ni^) ^« « 

first-order discretization scheme. 
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Achieving higher-order schemes for the approximation of L'^^-^ can be eas- 
ily obtained from the proof of Lemma 6.10. In fact, the required conditions 
are as foUows. Assume that for p>2, 

(6.15) / \y\H{yr\l + (|y|2/(y)-iA,t)(*'"2)/2)z.(dy) < Ct, 
J\y\<e 



(6.16) / \yVy{dy) < Ct\ 

Ay\<e 

For L'^_^^^, the idea used in the previous scheme is that the probabihty of 
having more than one jump in an interval of size T/n is of order (T/n)'^ and 
therefore they can be neglected if the goal is to achieve a scheme of order 
1. Obviously, in order to obtain a higher-order scheme, one has to consider 
the possibility of more jumps per interval. As an example, we consider the 
case of at most two jumps per interval with localization 1 = 1. 

For L^_^_l one can do the following: let Gg^dy) = C~^l|j^|>£i^(d2/), Cg = 
J\y\y£^{dy) and let Zf, ~ Ge independent between themselves and let Sf 
and Sf be two independent Bernoulli random variable independent of Zf, 
Z'k. Then consider the following cases. If 5f = define ^(x) = x, if Sf = 1 
and 5*2 = then '^(x) = x + h(x)Z\ and finally if 5f = 1 and S"! = 1 then 
Xl\x) = X + h{x)Zl + h{x + h(x)Zl)Zl. 

Define 

p, = P[5f = l](l + P[5| = l]), 

qg = P{Sl = \\P{S\ = \\. 
In this case we have the following lemma. 

Lemma 6.13. lfC-^P\S\ = 1, 5f = 0] < Ct and C-^P[Sf = 1, 5f = 1] < 
Ct then Assumption (M ) is satisfied with 

E[\x'4,ix)\P]<{l + Kt)\x\P + K't 

for all p>2. 

Proof. The result follows clearly from {f{x) = |x|^) 
P[St = 1, SI = 0]\E[f{x + h{x)Z') - fix)]\ 

<Ctil + \x\P)(l+ I \y\Pu{dy)\ 
P[Sl = 1, 5| = l\\E[f{x + h{x)Zl + h{x + h{x)Zl)Zl) - f{x)]\ 



<Ct{l + \x\'P)(l+( [ \y\'Pv{dy) 
V \J\y\>e 



□ 
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Lemma 6.14. Assume that \C~^Pe - t\ < Ct^ and \2C^'^qe - < Ct^ 
then 



E[fiXtn] - fix) - tLl^if{x) - -^{Li+,Yf{x) 

<Ct^\\f\\ci{l + \xr^){l+{j \y\u{dy) 
" V \J\y\>e 



Proof. As before let f e Cp then 

= fix) + / {fix + h{x)y) - f{x))u[dy)C-^P[Sl = 1, 5| = 0] 
J\y\>e 



+ E 



fix + hix)y + hix + - fix)vidy) 



\>e 



xC7ip[5f = 1,51 = 1] 
fix) + Ll^J{x)C~^P[Sl = 1, 5| = 0] 

fix + hix)y + hix + hix)y)yi) - 



+ 



|y|>e"'|2/|>£ 
2 i 



xC-'P[Sl = l,Sl = l\ 
= fix) + Ll^Jix)C~\P[St = 1] + P[St = 1, 5| = 1]) 
+ (Ll+i)2/(x)C-2p[5f = 1,51 = 1]. 
Therefore 



Eifixf'^)] - fix) - tL^Jix) - -iLl^.ffix] 



<|L3+i/(x)||C7-V-t| + l(ii+i)V(^) 
FinaUy note that 
{Ll+i)'fix) 

Ld+if{x + hix)y) - Ll^Jix)uidy) 



r-^n - — 



e<\y\ Js<\yi\ 



ifix + hix)y + hix + hix)y)yi) 

- 2fix + hix)y) + fix))vidyi)uidy) 
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/ / / y f {x + h{x)y + ah{x + h{x)y)yi)h{x + h{x)y)yi 

Je<\y\ Je<\yi\ JO 

— Vf{x + ah[x)y)h{x)y dav{dyi)u[dy) 

1 

V/(x + h{x)yi + ah{x + h{x)yi)y) 

X / Vh{x + (3h{x)yi)h{x)yidf3ydav{dyi)v{dy) 
Jo 



e<\y\ Js^yil Jo 



+ I I II D^f{x + ah{x)y 

Je<\y\ Je<\yi\ Jo Jo 

+ P{h{x)yi + a{h{x + h{x)yi) - h{x))y)) 

X h{x)yi 

+ a(^j Vhlx + ^h{x)yi) d'j h{x)yi^y,h{x)y dfi dav{dyi)v{dy). 
This finishes the proof. □ 



A similar statement can be achieved if we limit the number of jumps in 
any interval. The parallel of Theorem 6.6 can also be stated in this case. 

6.5. Example: tempered stable Levy measure. Now we consider the pre- 
vious approximation in the case that the Levy measure v defined on Rq is 
given by 

z.(dy) = ^i^(c+e-^+l^ll,>o + c_e-^-l^ll,<o)dy. 

The Levy process associated with no Brownian term and the above Levy 
measure v is called by: 

• Gamma: A+, c+ > 0, c_ = 0, a = 0. 

• Variance gamma: A+ , A_ , C-|_ , c_ > 0, a = 0. 

• Tempered stable: A+jA-, c+,c_ > 0, < a < 2. 

Then, we have that for a € [0, 1) 

/ |y|M'^y)~e^-", k>l. 
J\y\<e 

Then sup^g^g,!] i|y|<£ < co. For L^^;^, we consider as localization 

function l{y) = \yY , then the conditions of Lemma 6.10 are satisfied if a < 
r < 2 and e = ti/(3-a). 
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For L^^^, we consider as localization l{y) = 1, then Lemma 6.12 is sat- 
isfied, for example, in the following case. Let P[S^ = 1] = e^'^^'^l^.*) where 
Ce ~ a{e,t) = -e"log((t2 + as e = tV(3-a) ^^^^ ^ime that 

a = log((t + l)i(3-2a)/(3-a))_ 

In the case of Lemma 6.14, one choice of parameters is 

P[Sl = 1] = i(6-3a)/{3-a)(^ + 1)(1 + T/^^"")), 
^[^2 = 1] = 2(l + p/{3-a))- 

The choice of r in the above scheme is related with variance/importance 
sampling issues. 

Final comment: In this article we have presented a general setup to han- 
dle what maybe called operator decomposition methods. In particular, the 
method is useful when considering approximations of expectations of func- 
tional of diffusions (for another similar approach, see Alfonsi [1]). The ap- 
proximation problem is divided in components, each one driven by a sin- 
gle process. This single process, called the coordinate process can be ap- 
proximated to a high order using an appropriate (stochastic) Runge-Kutta 
scheme if the driving process is the Brownian motion. In the case that the 
driving process is a Levy process one can decompose the Levy measure in 
various pieces to facilitate the analysis. Note that sometimes is not needed to 
know how to simulate Y but only the functional form of the Levy measure. 
In comparison with the proposal presented in [9] , where high-order multiple 
integrals driven by different Wiener processes have to be simulated at each 
step, we believe that the present methodology is a better scheme. 

The issue that local approximations of high order are interesting to study 
in comparison with Romberg extrapolations as introduced in [19] is similar 
to the discussion of using Runge-Kutta approximations in comparison with 
Romberg extrapolations to approximate solutions of ordinary differential 
equations. We believe this article helps to open the path in this direction. 
In fact, it is somewhat clear from Theorem 4.2 that the leading constants in 
a Euler+Romberg method and a Runge-Kutta method do not coincide. 

Finally, we used the structure of this construction to easily introduce 
and analyze the asymptotic error of an approximating scheme for solutions 
of stochastic differential equations driven by Levy processes with possibly 
infinite activity. 

APPENDIX 

In this section we assume condition {A4p). 
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Lemma A.l. Let f = fs{x) £ C^'^{[0,T]xR^). Then a map s ^ Psfs{x) 
is Lipschitz continuous for all x S R^. 

Proof. Note that 

\Ptft{x) - Psfs{x)\ < \Ptft{x) - Ptfs{x)\ + \Ptfs{x) - Psfs{x)\. 

Using the Lipschitz properties of 1 1— > ft{x) and 1 1— > Ptfs{x), the proof fol- 
lows. □ 



Lemma A. 2. Let g -.[OjT] be a Lipschitz continuous function. Then 
we have 



(A.l) 



— 9{kT/n) - / g{s) ds 



< 



C{T,g) 



n 



Proof. From the assumption we immediately obtain 



2^ rkT/n 

—g{kT/n) - / gis)ds 



— 9 ' 



rkT/n 
l{k~l)T/r. 

where C depends on T and the Lipschitz coefficient of g. This implies (A.l) 
□ 



A.l. Some properties of Levy driven SDEs. We start with the differen- 
tiability properties of Xt{x) in x. The following material can be found in 
[7, 8, 10, 17] and [18]. We quote them here for completeness. 

Lemma A. 3. There exists a version of Xt{x) such that a map x i— > Xt{x) 
is infinite times continuous differentiable almost surely and in the -sense. 
Moreover, we have for p>2, 



E 



sup \Xt{x)\P 

0<t<T 



(A.2) 

and 
(A.3) 

for any multi-index a with \a\ > 1. 



<Cip,T){l + \x\P) 



sup E 



sup \d^Xt{x)\'P 

0<t<T 



< oo 



Proposition A. 4. Let f e with p>2. 
(i) Then Ptf ^C^ for aUt>0 and 



(A.4) 



sup ||PJ||c,"^<C||/||c7,. 
te[o,T] 
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(ii) Ifm>2, then Lf G Ci;!^2^ and 

\\Lf\\c„._-2<C\\f\\c^. 

(iii) /// G Ci'™([0,r] X R^), then {dtLf){t,x) = {Ldtf){t,x). 

Proof. The proof of (i) follows by interchange of derivation and ex- 
pectation together with the moment estimates in Lemma A. 3. Recall that 
L = Ef=o as defined in (2.5). 

(ii) We only do the proof for Ld+i- We have 

{fix + h{x)y) - fix) - Vf{x)h{x)T{y))v{dy) 

1 



< 



Vf{x)h{x){y-r{y))u{dy) 



+ 



llo de^f('^ + ^^('^)y)dOi^{dy) 



<c\\f\\cAi + \x 



□ 



Proposition A. 5. Let f G Cp. Then Pt and L are commutative and 
Uf{t,x) ■.= Ptf{x) is the solution of the integro- differential equation: 



d_ 

lit 



Uf{t, x) = Luf{t, x), 



Uf{0,x) = f{x). 

Let / G Cp'^^'^. Then the commutativity of Pt and L implies that L'^Uf 
(= ULmf) is differentiable in t and the solution to similar integro-differential 
equations. That is, 



d_ 
di 



{L"'uf){t,x) = L{L"'uf){t,x), 



(L™n;)(0,x) = (L™/)(x) 
for each m > 0. Consequently, applying Taylor's expansion to tij, we have: 



Proposition A. 6. For f e Cp™+^, 



Ptf{x) = Y^-L^f{x) + 



fc=0 



it - sY 

ml 



Ps{L"^+'f){x)ds. 



l,m-2 



Furthermore, if f ^ C™ with m>2. Then Ptf ^ 
Summarizing this section, we have 
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Corollary A. 7. Ptf{x) = E[f{Xt{x))] and Q\f{x) = E[f{Xi{x))] {i = 
0, satisfy conditions (M) and TZ{m,t^). That is, for p G'N , 

E[\Xtix)\'^P] < (1 + Kt)\x\^P + K't 

for some constant K = K{T,p), K' = K'{T,p) > and 

J<m{Pt) = ry^^' 

J<m{Q\) = 2Z Ti^i 

for any m G N. 
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